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Abstract 

We discuss the gamma-ray absorption in the inner region of the microquasar SS433. Our investigation includes several con- 
tributions to the opacity of this system. They result from the ambient fields generated by the primary star, possibly an A-type 
supergiant, and a very extended disk around the black hole. Besides the sharp and dramatic absorption effect that occurs every 
time the star crosses the emission zone, we find in the UV photon field from the extended disk an important source of absorption 
for very high energy gamma-rays. This results in periodic gamma-ray observational signatures. 
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1. INTRODUCTION 

SS433 is a particularly interesting X-ray binary system 
chat consists of a donor star feeding mass to a black hole [1] 
in orbit with a period of 13 d. From the vicinity of the com- 
pact object two oppositely directed jets are launched devel- 
oping regular precession with a period of 162 d. These jets 
can be considered as 'dark' [2] because their main power 
output is given by their kinetic luminosity, ~ 10 39 erg 
s _1 . Indeed, the ejected matter has determined the defor- 
mation of the nebula W50 that surrounds the SS433 system 
[3]. 

Gamma-ray emission from similar objects has been con- 
firmed very recently [4,5,6], so the study of absorption in 
the complex case of SS433 is important to characterize the 
possible high energy emission that this source might present 
to instruments such as GLAST and the new Cherenkov ar- 
rays MAGIC II and VERITAS. For a recent discussion of 
the absorption in the case of other systems like LS I +61 
303 and LS 5039, see for instance, Ref. [7,8,9]. 
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Regardless of any specific gamma-ray emission process 
that could operate in SS433, in the present work, we shall 
focus on the possible sources of gamma-ray absorption tak- 
ing place mainly in the inner region and in the immediate 
vicinity of the binary. There, the interactions with photons 
and matter can be important and cause significant absorp- 
tion. In particular, if the emission is originated near the 
compact object, we find a dramatic absorption effect oc- 
curring every time the companion star eclipses the emis- 
sion zone. Furthermore, at very high energies (i? 7 > 50 
GeV) , we predict a very important absorption effect caused 
mainly by the UV emission from the extended disk around 
the accretion disk, and secondly by the mid-IR emission 
from the same disk. 

We have based our work on the set of parameters that are 
currently believed to describe the system SS433 after more 
that 20 years of observation and study. These parameters 
are summarized in the next section. The resulting absorp- 
tion signatures are presented in Sect. 3, and a discussion is 
left for Sect. 4, where we analyze the possible detection of 
the specific absorption features. 
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2. THE SOURCE 



The mass loss rate in the jets of SS433 is mj = 5 x 



3. OPACITY TO GAMMA-RAY 
PROPAGATION 



10 7 M yr 1 and their bulk velocity is « 0.26c. The Attenuation of the putative gamma radiation emitted by 



normal to the orbital plane makes an angle k 21° with 
the approaching jet and an angle i = 78° with the line of 
sight. The line of sight then makes a angle ij with the jet 
which is time-dependent because of precession (see Fig.l). 

A thick expanding disk wind believed to be fed by the 
supercritical accretion encloses also the star [10]. According 
to Rcf. [1] , the disk wind has a half opening angle a w ~ 30°, 
a mass loss rate M w ~ 10~ 4 M Q yr _1 , and a velocity v w <~ 
1500 km s _1 . Emission of UV photons from this extended 
disk will be considered using a blackbody distribution with 
T uv = 21000 K following Ref. [11], while mid-IR emission 
will be taken into account using a fit for the flux density 
presented in Ref. [12]. 

The spectral identification of the primary star has been 
difficult because of the presence of the extended disk, since 
the star is often partially or totally obscured by it. After 
convenient observations at specific configurations of pre- 
cessional and orbital phases it has become quite clear that 
the star is an A-supergiant [13,14,15]. We assume the mass 
of the components as derived from INTEGRAL observa- 
tions [13], Mbh = 9M and = 30M© for the black hole 
and the star respectively, which correspond to an orbital 
separation a ~ 79 Rq for a zero-eccentricity orbit as is the 
case for SS433. Since the star is believed to fill its Roche 
lobe, the implied radius according to Ref. [16] is = 
0.49g- 2 / 3 a[0.6g- 2 / 3 +ln(l + g - 1 /3)]- 1 ~ 38jRq; where 

q = Mbh/M*. We shall consider the emission of soft pho- 
tons from the star through an associated blackbody distri- 
bution with an effective temperature T* = 8500 K [14,15]. 
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Fig. 1. Schematic view of SS433. The approaching jet is most of the 
time closer to the line of sight and the receding one is oppositely 
directed. 



this source would be mainly due to e~e + pair production 
arising from 77 annihilation, particularly with lower energy 
photons coming from the extended disk and from the star. 
Matter, on the other hand, also offers a significant target 
for gamma-rays emitted by the system. We shall calculate 
also the effect of 7 TV interactions around the star and the 
extended disk. 

When a gamma-ray of energy E 1 travels a distance dp 1 
in a photon field, there is a differential optical depth asso- 
ciated with it. Absorption correspondingly arises from the 
gamma-ray interaction with soft photons of energy E. As- 
suming that the former runs following e 7 and the latter are 
directed along e pn = (sin 9' cos </>', sin 8' sin </>', cos 8'), the 
corresponding optical depth is given by [17], 



dr 77 = (1 — e 7 ■ e p h)n p hcr 77 dp 1 dE dcosd' dej) 1 . 



(1) 



Here, n p h is the density of the soft photons per solid angle 
and soft photon energy units. The angles 9' and <j)' are taken 
in a convenient coordinate system and the cross section for 
the process 77 — > e + e~ is given by 



2£(£ 2 -2) + (3-£ 4 )ln 



1+j 



where ro is the classical electron radius and 
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3.1. Optical depth due to the companion starlight 

Starlight photons coming from the companion mid-A su- 
pergiant at a temperature T ~ 8500 K [14,15] can cause 
absorption of gamma-rays. To see this clearly let us con- 
sider some simple geometry. 

In Fig. 1 the observer is assumed to lay in the xz-plane 
and the gamma-ray path is described by the vector 



Pi 



p 7 (xsini + zcosi) = p 7 e n 



(4) 



with i = 78° as mentioned. The position of the star is given 
by a = a(x cos 4>+y sin (f) and we suppose that the gamma- 
ray is produced in the jet at 



Zj = Zj (x sin 9 cos ip + y sin 9 sin tp + z cos 8) , 



(5) 



where ip is the precessional phase. The position where the 
interaction with the soft photon takes place is indicated by 



Try Zj ~\~ . 



(6) 
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In this case, we set ^ p h = (r 7 — a) as the vector connecting 
the center of the star with the gamma-ray, and we call a 
the angle that it makes with e 7 . 

At each particular gamma-ray position, a convenient ref- 
erence system can be chosen to perform the integration 
in the angles 9' and <j>' in expression (1). This is done, as 
in Ref. [18], orientating the z'-axis along the vector Z p h- 
Accordingly, the angle <j>' varies between min — and 
'/'max = 27r and 9' varies between 9' min = and max = 
arccos(yT — (i?*/^) 2 )- Further orientating the rr'-axis so 
that the gamma-ray path is contained in the x'z'-plane, the 
gamma-ray direction can be written as e 7 = (sin a, 0, cos a) 
and hence, e 7 • e pri = sin a sin 9' cos </>' + cos a cos 9' . How- 
ever, the dot-product involved for the evaluation of the 
cross section (2) and the minimum photon energy (10) is ap- 
proximated by cos a, which does not depend on the primed 
angles. This last simplification that permits the analytic in- 
tegration on <j)' and 9' , means only that in order to compute 
the cross section and the minimum energy, the starlight 
photons are all considered as being originated at the center 
of the star. 

Integration in 9' and </>' of expression (1) thus gives 
dr 

— = 2Kn i ,(E)a rr {E,R r ) 

x [cos0 min -cos0 max (7) 

COS Ot - . o f.1 . 9 n/ \ 

— (sin 6> max -sin min ) . 

The density of radiation from the star can be approximated 
as usual by 



2E 2 



(h cf{e E / kT * - 1) 



(ph cm 3 crg x sr 1 ). 



(8) 



Further integrating in target photon energy E and along 
the gamma-ray path p 7 yields the resulting starlight con- 
tribution to the optical depth 



T* = 



where 



E m i n — 



dE I dp. 



dr+ 



dEdp~j 
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(10) 



Once one relates the position vectors with time, one can 
obtain t*(.E 7 , t) as shown in Fig. 2. These gamma-rays were 
assumed to be originated mostly near an injection point at 
the base of the approaching jet, zq = Rq/ ^an£ ~ 1.3 x 
10 9 cm, where ^ ~ 0.6 is the half opening angle of the jet 
(see Ref. [19]), and the initial jet radius is taken to be i?o = 
10GM bh /c 2 . It can be seen from Fig. 2 that this component 
of the opacity is notoriously more important when the star 
is in the foreground of the emission point, namely when it 
crosses the x-axis. 

We note that to perform the intgral in Eq.(9), we have 
assumed for simplicity that the radiation density given by 



Eq.(8) was artificially valid inside the star. However, as will 
be shown below, the absorption caused by 7./V interactions 
is by far the dominating one if the gamma-ray path directed 
to Earth has to travel through the star. This implies that, 
although the peaks in Fig. 2 could look different, an accu- 
rate expression for the radiation density inside the star will 
not affect the total optical depth. 




Fig. 2. Starlight optical depth. 



3.2. Optical depth due to the emission from the extended 
disk 

The extended disk wind is believed to be the origin of 
both the observed mid-IR and UV emission. In the first 
case, the reported free-free emission was detected in the 
range of wavelengths 2-12 /im for which we adopt the fol- 
lowing flux density fit [12], 



F ph = 2.3 x 10~ 23 ( — 
\/im 



-0.6 



(erg S - 1 cm- 2 Hz- 1 ), (11) 



which corresponds an emitting region of radius i? out = 
50i? Q . We need to estimate the radiation density of these 
mid-IR photons, i.e., the number of photons per unit vol- 
ume, per unit energy per unit solid angle spanned by the 
emitting region. The corresponding mid-IR photon flux per 
unit frequency arriving at Earth can be written as 



dN, 



P h 



dt dA dv 



%(ph cm-V 1 ^- 1 ) 
E 



(12) 



Considering that dE/dv = h, dl p h/dt — c, and that 
dlphdA = dV, the radiation density coming from the whole 
emitting region is 



dN, 



P h 



dVdE E he 



h (ph cm 3 erg 1 ). 



(13) 



Since the adopted fit for F p h corresponds to the particular 
precessional phase ^>ir = 0.31 x 2tt (see Ref. [12]), in that 



case, our line of sight makes an angle ijjR with the normal to 
the disk midplane such that cosijjR = sin i sin 9 cos V'ir + 
cos i cos 9. Then, the area of the emitting region 7ri?Q Ut spans 
a solid angle Aft — (it Rl ut cos ijjR)/d 2 as viewed from 
Earth, and we can estimate the corresponding radiation 
density as 



1 f dN ph \ 
Um ~ Afi \dVdEj 



E he 71-cosijjR. V-Rout 



d 



(14) 



(ph cm 3 erg 1 sr 1 ). (15) 



We proceed to integrate the mid-IR contribution of the 
optical depth tir, using the method described in Appendix 
A for a disk-like emitting zone. The inner radius of this zone 
is taken as the inner radius of the extended disk, i?; n ss 
2GMbh/«2>, as suggested in Ref. [20]. The obtained result 
is shown in Fig. 3 for gamma-rays originated at the base of 
the approaching jet. 
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Fig. 3. Mid-IR optical depth as a function of time and gamma-ray 
energy for a gamma-ray produced at a height zq on the approaching 
jet. 

As for the UV emission, it was detected in the range 
of wavelengths - (1000 - 10000)1 [11] and, as a first ap- 
proach, we consider an associated blackbody temperature 
Tuv = 21000 K for the disk in an edge-on state, with a 
corresponding radius of the emitting zone R out — 33i? Q . 
Therefore, we consider a radiation density for the UV emis- 
sion as the one given in expression (8) with T+ — > Tuv- 

The integration of this contribution to the total optical 
depth is also performed as it is described in Appendix A. 
The obtained contribution is shown in Fig. 4 as a function 
of time and gamma-ray energy. 

It can be noted from this plot and from Fig. 3 that the 
maximum absorption occurs near the two times t ~ 55 d 
and t ~ 110 d. Since the mid-IR and UV emissions are 
treated in this first analysis as being originated in the mid- 
plane of the extended disk, these times correspond to the 



cases when the gamma-rays originated at z$ and directed 
to Earth run almost parallel to the mentioned emitting 
zone. This implies that 77 collisions will occur at larger an- 
gles resulting in higher cross sections and smaller threshold 
energies along a significant part of the gamma-ray path. 
At intermediate times, 55 d < t < 110 d, the mentioned 
gamma-rays have to pass through the midplane of the ex- 
tended disk, and once they have done so, the collision an- 
gles are smaller, the threshold energy is larger, and hence 
the absorption is lower. 

As for the dependence with the gamma-ray energy, it 
can be seen from Figs. 3 and 4 that the UV optical depth 
becomes important for E 1 > 50 GeV, while the mid-IR 
contribution will be significant for E 1 > 500 GeV. 
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Fig. 4. UV optical depth as a function of time and gamma-ray energy 
for a gamma-ray produced at a height zq on the approaching jet. 



3.3. Optical depth due to jN interactions. 

The nucleons of the star and the thick extended disk 
can also absorb gamma-rays and photo-produce pions. We 
consider the cross section as in Ref. [21], 



0~ 1 N 



340 Aib for 200 MeV < E 1 < 500 MeV 



120 /ib 



for E^ > 500 MeV 



(16) 



where the first case corresponds to the single pion channel 
and the second case to the multi-pion channel. 

Considering the unit vector Z = (sin 9 cos ip, sin 9 sin tp, cos ( 
which is perpendicular to the midplane of the disk wind, 
and the angle it makes with the gamma-ray position, 
9z = arccos(Z • r 7 /r 7 ), the density of the extended disk 
wind at the point r 7 is estimated to be 



Pw(r^,9 z ) = 



e 



flz-(|-«w)]e[(| + ow)-(9 



Here, is the Heaviside step function and the solid angle 
element is related to the disk wind half opening angle a w 
by A£l = 47rsina w . 

As for the star, we suppose that it has a density given by 



P*(r) 



e(r-K), 



(18) 



where r represents the distance to the center of the star. 

The optical depth due to 7 A interactions with the men- 
tioned nuclcons can be estimated as 



T-yN 



CO 

'■}) = / ap 7 a jN - 



(P* + Pw) 



(19) 



The obtained result is shown in Fig. 5 as a function of 
time for 200 MeV < E y < 500 MeV and £ 7 > 500 MeV. 
As it can be seen from this plot, periodic peaks of high 
absorption appear when the star is in the foreground. We 
note that the absorption caused by the disk wind is rather 
low. In particular, it vanishes when the approaching jet is 
most open to us (for t < 12 d and t > 148 d ) because in 
those cases the gamma-ray path originated at the base of 
the jet never travels through the disk wind. 
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Fig. 5. ~fN optical depth as a function of time for 
200 MeV < £ 7 < 500 MeV (solid line) and E 1 > 500 MeV (dashed 
line). The sharp peaks arc due to the orbital modulation. 



3.4. Total optical depth 



Analyzing the different contributions, it can be noted 
that below E 1 = 20 GeV, the only relevant source of ab- 
sorption is due to the 7 A interactions with the nucleons of 
the extended disk and the star. 

If we add up all the contributions calculated above for 
the case of gamma-rays originated at a height z in the ap- 
proaching jet, we can obtain the total optical depth which 
is shown in Fig. 6 as a function of time and gamma-ray en- 
ergy. It can be seen clearly from this plot that the absorp- 
tion caused by 7A interactions with the nucleons of the 



star is important along all the range of gamma-ray energies 
studied (200 MeV - 20 TeV). 
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Fig. 6. Total optical depth for gamma-rays with as a function of time 
for E 7 = 200 MeV (solid line), E-, = 800 GeV (dashed line), E 1 = 5 
TeV (dotted line), and E 1 = 10 TeV (gray solid line). 



3.5. Dependence of the absorption on the gamma-ray 
injection point 

We shall now consider the cases where most of the 
gamma-rays are produced at a given height zj > zo along 
the approaching jet. The calculation of the various compo- 
nents of the optical depth can be performed as described 
in the previous section substituting Zo by zj . 

We show the obtained total optical depth as a function 
of Zj and the gamma-ray energy in Figs. 7 and 8 for t ~ 
162 d and t <~ 81 d respectively, representing the cases 
when the absorption is maximum and minimum. Since at 
these times the star is not in the foreground, the optical 
depth is basically the one corresponding to the UV and 
mid-IR emission from the disk wind. At t = 162 d, when 
the approaching jet is mostly open to our direction, the 
optical depth decreases with zj because as zj increases, the 
distance to be traveled by the gamma-ray directed to Earth 
over which 77 interactions are important decreases. 

At t — 81 d there is an initial rise of the optical depth with 
zj. This is because at this particular time, the approaching 
jet as was referred here for simplicity, is actually not ap- 
proaching the Earth but directed away from it. Then, as zj 
increases, the distance along which 77 collisions are rele- 
vant increases for the gamma-rays coming in our direction. 
However, for still higher values, zj > 10 11 cm, the optical 
depth decreases with zj since the density of UV radiation 
from the extended disk wind is lower. 

4. SIGNATURES ON A GAMMA-RAY SIGNAL 

In the present work we have not specified any particular 
gamma-ray emission process. Regardless of this, one possi- 




Fig. 7. Total optical depth at t = 162 d for gamma-rays originated 
with energy E 1 at a distance Zj from the compact object along the 
approaching jet. 




Fig. 8. Total optical depth at t = 81 d for gamma-rays originated 
with energy E 1 at a distance Zj from the compact object along the 
approaching jet. 

bility is that since the particle densities are expected to be 
highest in the inner jet, most of the gamma-rays could be 
produced at a short distance Zj from the compact object. 
In view of the recent detection of gamma-rays from simi- 
lar systems, such as LS I +61 303, LS 5039, and Cyg X-l 
[4,5,6], and in an attempt to give an idea of the effect that 
the obtained optical depth may cause on an out-coming 
gamma-ray flux, we can assume that most of the radiation 



is produced near the base of the jets (at zj ~ z ). In these 
conditions, if the spectrum of produced gamma-rays fol- 
lows a power law as 

J 7 = K^E- 2 (ph erg^sr^cm-V 1 ), (20) 

where K 1 is a constant, the total luminosity within an en- 
ergy range (£™ in , £™ ax ) is 

L 7 = AA J dE J E J J 7 (E 1 ), (21) 

where AA is the element of area of the emitting region near 
the injection point. It follows, then, that 

K 7 AA = — ^ssr, (22) 
In ' ■ 

1 

which allows us to estimate the corresponding photon flux 
to be detected at the Earth as 

$ ^ ) = 4^ / dE^EJe-^^). (23) 

B m.n 

For illustration, we assume that the equivalent isotropic 
gamma-ray luminosity between E™ m = 200 MeV and 
fi™ = 20 TeV is L 7 w 10 36 erg s~\ The resulting flux is 
shown in the upper panel of Fig. 9, whereas the flux corre- 
sponding to energies E 1 > 800 GeV is shown in the lower 
panel as compared to the upper limit given by HEGRA 
[22] . In the first case the obtained mean photon flux is 

( $ 7)U 7 >200MeV = 5 x 10" 8 P h cnTV 1 , (24) 
and in the second case, 

<*7>k> 8 ooGcV = 4 -7 x 10- 14 ph cm-V\ (25) 

which is below the HEGRA cut, $ij m = 8.9xl0" 13 ph cm _2 s _1 . 

We emphasize that these values are assumed only to 
make a qualitative description of a possible gamma-ray sig- 
nal and that their true values may be obtained in a more 
detailed study regarding the emission process that could 
operate. 

5. DISCUSSION 

As it is reasonable to expect, the most noticeable ab- 
sorption signature that can be imprinted on a gamma-ray 
flux from the inner regions of SS433 is given by the reg- 
ular and dramatic absorption caused by the star when it 
crosses the line of sight. This effect takes place approxi- 
mately once every 13 days during the time that the star 
blocks the emitting region (« 2d). In practice, this could 
serve to determine the size of the companion star as well 
as the mass distribution of the stellar atmosphere as long 
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Fig. 9. Gamma-ray flux a function of time for Ej > 200 MeV (upper 
panel) and for E-y > 800 GeV (lower panel) as compared to the upper 
limit given by HEGRA. The assumed gamma-ray luminosity emitted 
is Lj = 10 36 erg s _1 and the maximum energy is B™" = 20 TeV 

as sufficient time resolution can be achieved in the 7-ray 
lightcurves. It seems possible that, given the little absorp- 
tion that corresponds to energies below ~ 50 GeV and the 
fact that more events are expected for lower energies, the 
GLAST instrument could detect a signal possibly with the 
time behavior predicted here. Another detectable absorp- 
tion feature to observe at these energies corresponds to the 
7iV interactions with the nucleons of the extended disk. 
This is expected to cause, because of the precession, a long 
term modulation of the signal reducing the maximum flux 
by a significant fraction as shown in the upper panel of Fig. 
9, provided that the bulk of the emission is produced near 

Zq. 

As for higher energies, absorption through 77 interac- 
tions also becomes important. The dominant contribution 
is caused by the UV emission from the extended disk for 
^7 It 50 GeV, while the mid-IR photons there originated 
can cause rather significant attenuation for E 1 > 500 GeV. 
The absorption due to the starlight photons is significant 
when the star is nearly in the foreground, almost eclipsing 
the black hole. This can be seen in the lower panel of Fig. 9, 
where the gamma-ray signal decreases and increases quite 
smoothly before and after the star has actually eclipsed the 
emission zone. 

The fact that HEGRA has failed to detect anything 
greater than $ij m = 8.9 x 10~ 13 ph s _1 cm" 2 for E 1 > 800 
GeV, implies that maybe SS433 is not intrinsically lumi- 
nous at these energies, or that gamma-rays could be actu- 
ally produced at an inner injection point as illustrated here 
with a given luminosity, but nothing could be detected at 
the Earth because of absorption. If the assumed luminos- 
ity in gamma-rays is similar to the real intrinsic luminos- 
ity, a detection at these energies seems more difficult but 
not impossible with the forthcoming Cherenkov telescopes 
MAGIC II and VERITAS. 

Finally we remark again that the conclusions drawn in 
the present study are independent of the gamma-ray pro- 



ducing mechanism. Taking into account reprocessing of 
gamma-rays due to cascading effects [23], would not ci- 
ther alter the absorption patterns here predicted. In case of 
an hadronic mechanism for gamma-ray production in the 
source (e.g. [24]), neutrino emission should be expected, as 
pointed out by several authors [25,26]. Contrary to gamma- 
rays, the neutrino signal is not affected by absorption, and 
its detection by instruments like IccCube might yield light 
on the relativistic particle content of the inner jets. 

A comprehensive study, including the description of the 
physical mechanisms allowing the high-energy emission 
processes in SS433 will be presented elsewhere. 
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Appendix A. Calculation of the gamma-ray 
absorption near a disk-like emitting region 

We consider the absorption of gamma-rays as they travel 
near a disk that emits soft photons. This disk, with an 
inner radius R m and an outer radius i? ou t, lies in a plane 
perpendicular to the jets, where the gamma-rays emerge. 
For gamma-rays produced with energy E 1 , traveling in a 
direction e 7 , the differential optical depth due to photons 
with a radiation density n pn and a direction given by e p h, 
can be written as (e.g [27]) 



C?T 77 — (1 6ph ' 6 7 )?7< p hC^ 77 



cos r)RddRdd<f)d 



dp-ydE. (A.l) 



Here the Rd is the length of the vector Rd signaling a point 
in the disk, <pd is its corresponding azimuthal angle, and 
Z p h is the magnitude of the vector l p h connecting the point 
on the disk to the gamma-ray position. This vector makes 
an angle j) with the normal of the disk plane. All the angles 
and vectors involved in (A.l) are then to be expressed in a 
new coordinate system with its Z-axis oriented along the 
axis of the approaching jet. 

It is worth noting that the radiation density n ph appear- 
ing in expression (A.l) represents the number of soft pho- 
tons per unit volume, per unit energy, per unit solid an- 
gle. Hence, for photons originated on a differential area el- 
ement da = RddRddfid on the disk plane, the differential 
solid angle spanned is dQ = (da cos i]) / l^ h , which appears 
in expression (A.l). 

In the system fixed to the compact object the unit vector 
e 7 is given by Eq. (4) in terms of the unit vectors x , y 
and z, and the basis of the new system can be obtained by 
transforming them according to 



( cos ip cos 9 — sin ip cos ip sin 9 ^ 



R 



sin ip cos 9 cos ip sin ip sin 6 ■ 
— sin 9 cos 9 



(A.2) 



Therefore, the unit vectors of the new coordinate system 
can be expressed in terms of fixed basis as 



X = x cos 9 cos ip + y cos 9 sin ip — z sin 9 

Y = —x sin ip + y cos ip 

Z = xsin9 cos ip + y sin 9 sin ip + z cos 9, 

and the unit vector <=L can be written as 



(A.3) 
(A.4) 
(A.5) 



Y e-v • Y 



Z ( ely • Z 



e 7 = X (cos 9 sin ip sin i — sin 9 cos i) 

— Y (sin -0 sin i) (A. 6) 

+ Z (sin 9 cos ip sin i + cos 9 cos i) . 
Since the position of the gamma-ray is 

r 7 = e 7( o 7 + Zzj 

and the position on the disk is 

Rd = Rd(xcos(pd + Ysiri(pd] , (A.7) 



the vector / p h can be obtained as l p h = r 7 — Rd- 

The resulting optical depth is given by the quadruple 
integral 



2tt R a 



J dp 7 J d(p d J dR d J dE 

_R in £ min 



dr. 
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dp^dcpddRddE 



where 



E, 



2(m e c 2 ) 2 



(A.8) 



(A.9) 



Ej(l - e 7 • e ph ) 

The integral can be performed using a Monte Carlo 
method, that is, introducing the variables x p , x^, xr, and 

Xe as 



p 7 = pix p 
<pd = 2nx<p 

Rd = Rin + (Rout — Rin)xR 

E = E m [ n -(- (E max E m \Yi)x e i 
then the integral (A.8) can be written as 



(A.10) 
(A.11) 
(A.12) 
(A.13) 



dx p \ dx^ / dx R / dx E f(xp 1 x ( p,x R ,x E ), (A. 14) 



o 



where 



/(xp^x^^x^xe) = 2irp 1 (R ou t - R in ) x 



(Emax E m [ n ^ 



dr-y 



dp~/d(pddRddE 



(A.15) 



Here the upper limit of the integration in p 7 is taken as 
pi ~ 10i?out where the integrand is significative. 

The method relies on the fact that, by the Mean- Value 
Theorem, the integral can be approximated by 



1 N 



(A.16) 



where each variable takes a random number between and 
1. The right hand side is then the statistical average of N 
evaluations of the function /. The error goes like l/\/N, 
which in this case of four integration variables makes it a 
quite accurate method as compared to the iterative ones 
(see e.g. [28]). 



